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The one-loop anomalies of a (i-dimensional quantum field theory can be computed by 
evaluating the trace of the path integral Jacobian matrix J, regulated by an operator 
exp(— /37?.) and taking the limit (3 to zero In 1983, Alvarez-Gaume and Witten 
made the observation that one can simplify this evaluation by replacing the operators which 
appear in J and TZ by quantum mechanical operators with the same (anti) commutation 
relations. By rewriting this quantum mechanical trace as a path integral with periodic 
boundary conditions at time t = and t = (3 for a one-dimensional supersymmetric non- 
linear sigma model, they obtained the chiral anomalies for spin ^ and | fields and self-dual 
antisymmetric tensors in d dimensions. 

Some time ago it occurred to us that one can also apply these ideas to the trace anom- 
alies. In a recent paper the first author proposed a bosonic configuration space path 
integral for a particle moving in curved space, and found the corresponding Hamiltonian 
71 from the Schrodinger equation. He exponentiated the factors ^/g in the path integral 
measure by using scalar ghosts, and obtained the trace anomaly for a scalar field in an 
external gravitational field in d = 2. In this article, we treat the general case of trace 
anomalies for external gravitational and Yang-Mills fields. We do not introduce a super- 
symmetric sigma model, but keep the original Dirac matrices and internal symmetry 
generators in the path integral. As a result, we get a matrix-valued action S. Gauge 
covariance of the path integral then requires a definition of the exponential of the action 
by time-ordering. The computations are simplified by using Riemann normal coordinates. 
We also replace the scalar ghosts by vector ghosts in order to exhibit the cancellation of 
all divergences at finite f3 more clearly. Finally we compute the trace anomalies in d = 2 
and d = 4. 
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1. Introduction 

The one-loop anomalies in local symmetries of a d-dimensional quantum field theory 
can be computed, as first demonstrated by Fujikawa by evaluating the trace of the 
Jacobian matrix J that arises when one varies the integration variables of the path integral 
under this local symmetry. This trace has to be regulated because the Jacobian is an infinite 
dimensional matrix, and as regulator one may take an operator of the form exp(— The 
anomaly is then obtained by letting (3 tend to zero. In general, there will be divergent terms 
proportional to powers of and a finite remainder. The latter is the anomaly, while the 
former are discarded by using several regulators whose /3's satisfy suitable relations as in 
Pauli-Villars regularization. The justification for discarding the divergences in anomalies 
is that the effective action is made finite by adding suitable counterterms so that its 
variation should also be finite. (If one regulates the quantum field theory with dimensional 
regularization, the counterterms which make the one-loop effective action finite produce 
also anomalies [Q, which are finite since the one- loop pole is cancelled by a zero due to 
classical gauge invariance). 

The choice of a regulator TZ is important. For some cases one knows from computations 
a regulator which yields consistent anomalies, or a regulator which yields a covariant 
anomaly, but only for consistent anomalies has a general theory been constructed 0. 
This theory provides a method to construct regulators which: i) preserve certain local 
symmetries at the quantum level, and ii) produce anomalies in other symmetries which 
satisfy the consistency conditions. For Einstein symmetries (general coordinate invariance) , 
Fujikawa already determined the regulators which satisfy (z). For covariant anomalies, no 
general theory is known as yet, although in practice one usually can guess what these 
covariant regulators will be. In our work below, we are in the happy circumstance that 
the consistent regulators which preserve the Einstein symmetry are actually covariant, as 
we shall discuss. 

The local symmetries we shall consider are local scale symmetries. The corresponding 
anomalies are called Weyl anomalies, or trace anomalies because they arise as a non- 
vanishing trace of the stress tensor T^j^ at the quantum level. For example, consider a real 
scalar field ^, which couples to gravity in d dimensions with action 

(1.1) 

where ^ = 4^1^ • We shall always work in Euclidean space so that g = det g^jj. This 
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S = j dx^ ^ 



classical action is invariant under the following Weyl rescalings 



2-d 



g^iu g^iu' = ^9fiu; cj)^ ^' = vt 4 (1.2) 

The general theory of consistent regulators states that a consistent regulator without Ein- 
stein anomalies is obtained by using (f) = (f) a,s quantum field. Clearly 

Sw(i>=^(^(t>: (1-3) 



where a = logO is taken to be infinitesimal. Rewriting the action in (|1 . 1| ) in term of (j), 
the following regulator is obtained 

7^™-^ = -g'-^d^g-^g^^'d^g--^ - ^R. (1.4) 

This is clearly not an Einstein scalar. Under a rescaling of the external metric in the path 
integral Z[g] = J V(f)exp{—h~^S[g,(f)]), followed by a compensating rescaling of (f) as in 
( |T73|) such that S remains invariant, the following regulated Jacobian is produced: 

[ rfxV^a(x)^7'^^(x)(T^.) = 2nlimf^exp(-/37^^°'^^)Y (1.5) 

where we have defined the stress tensor by 

Since the Jacobian contains no derivatives, we can use the cyclicity of the regulated trace 
to cycle the factor g~'^ on the right-hand side of 'JZ^°'^^ in (|1.4| ) to the left-hand side, 
obtaining the trace anomaly of a scalar field in d dimensions 

Alf^ = I dx^ <y{x){T) = niim Tr a(x) exp(-/37^[)°") 



where the subscript zero indicates that we are dealing with spin zero fields and where T 
denotes the trace of the stress tensor. For Einstein anomalies, one cannot replace TZ^"^^ 
by 'R.'^"'" as the Jacobian contains derivatives |^ . The reason that this is possible for Weyl 
symmetries follows already from the well-known fact that the consistent Weyl anomaly is 
Einstein-covariant . 
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The problem is now to evaluate this trace. Fujikawa and others have used 
a complete set of plane-waves and integrated over plane-wave momenta, but the expres- 
sions become extremely unwieldy. DeWitt did develop, much earlier, a heat kernel 
analysis, according to which an ansatz for the answer is made, which is then verified by 
explicit computation order by order in (3. In the ansatz, "bilocal" tensors appeared (such 
as the Synge world function a{x,y), which measures the square of the distance between 
two points) and a calculus for the left and/or right derivatives of such objects was needed. 
Of course this approach yields much more than only the anomalies, for example one ob- 
tains in this way the complete one-loop effective action. Yet, a simpler method, which 
may only give the anomalies, might be welcome. Some time ago, Alvarez-Gaume and 
Witten PI observed that one can in general compute anomalies in a much simpler way by 
replacing the d-dimensional operators in J and TZ, such as x^, the Dirac matrices 
7^^, the metric g^jj{x), the connection A^{x) and the corresponding generator of internal 
symmetries {T°-)i^ , by quantum mechanical objects q^(t) and p^{t), 4'^{t), g^vio), ^m(?) 
and c*'^{t){T°')i^ Cj{t), where c** and Cj are anticommuting "ghosts" which contract the 
indices of internal symmetry generators such as Yang-Mills or Lorentz matrices. All that 
is needed is that they satisfy the same anticommutation rules. They then considered chiral 
anomalies and rewrote the trace of the quantum mechanical operator Jexp(— /37^), with 
J = 7^ as a one-dimensional path integral 

Z=l Vq^Vt/j^VciVc*' exp( - dt C^^'^X (1.8) 

^^^S 1°' '-t \ Jo J 



ABC for c,c* 



where C^'^^ is the Lagrangian for a one- dimensional supersymmetric non-linear sigma 
model. The (i-dimensional chirality operator 7^ became the fermion number operator 
(— )^, because (— )^ anticommutes with -0^ just as 7^ with 7^. Since (— )'^ commutes with 
and the ghosts, the tj}^ fermions acquired periodic boundary conditions (PBC), whereas 
the c** and Cj kept their antiperiodic boundary conditions (ABC). A simplification was due 
to the topological nature of the chiral anomaly: the path integral itself is /3-independent 
and hence could be evaluated by taking /3 small, and keeping only non-vanishing terms. 
In particular, no divergences proportional to (3~^ could occur. For the actual evaluation, 
Alvarez-Gaume and Witten expanded q^{t) and '(/''^(t) in fluctuations about a classical 
constant solution q^(t) = q^ and ijj^it) = i/'q , and only the terms quadratic in fluctuations 
did contribute. The result for the anomaly could be written in a covariant form by using 
normal coordinates and, since they used a covariant regulator, this covariantly looking 
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form was also valid in general coordinates. 

In this article we shall consider the trace anomalies, and extend the work of Alvarez- 
Gaume and Witten to this case. There are, however, important differences: 

i) The trace anomaly has no topological meaning, and as a consequence its path 
integral representation does depend on the regulating parameter (3. We must therefore 
expand in terms of /3, which amounts to a loop expansion on the world line with coordinate 
t. For the one-loop trace anomaly in d dimensions we shall need f + 1 loops on the world 
line. Hence the algebra becomes more complicated than for chiral anomalies and subtle 
problems which did not need to be addressed in p|, now have to be solved. For example, 
the naive relation between the Hamiltonian (=regulator) and the corresponding Euclidean 
action in the path-integral breaks down, and extra terms proportional to the curvature 
appear. The occurrence of such terms was first noticed by DeWitt 0. (In fact, in his well- 
known 1948 article 0, Feynman had already indicated that for models such as non-linear 
sigma models, it might not be correct to use the classical action and classical paths of a 
free particle to define the path integral.) 

ii) For the case of spin ^ fields, our Jacobian has no factor 7^, hence no operator (— )^; 
therefore all fermions will now acquire antiperiodic boundary conditions. Hence, there are 
no constant modes i/^q about which to expand. This, in itself, is no problem: we could 
expand ifj^ only into antiperiodic fiuctuations. 

iii) As pointed out in [||, the anticommuting creation and absorption operators c** 
and Ci should only act on the one-particle states in the quantum-mechanical trace of the 
operator Jexp(— /37^) in order that c*^T°'i^Cj acts in the quantum-mechanical problem in 
the same way as the matrices did in the original problem. Thus we would need a 
quantum-mechanical projection operator, which should become an extra internal kernel in 
the corresponding one-dimensional path integral. Although we have found this kernel, see 
(B.8), the boundary conditions to which it corresponds lead to complicated propagators, 
and therefore we decided to follow a different path. We keep at all times the Dirac matrices 
7'^ and internal symmetry generators T" as matrices, and never introduce any i(;'^(t), c**(t) 
or Ci{t) at all. Although to our regret supersymmetry leaves the stage by this approach (it 
could have been present but for the spin ^ regulator, and it would anyhow have been broken 
by the boundary conditions), it seems the simplest way to compute the path integral. In 
fact, we shall rewrite the action S in terms of a connection and a potential V, which 
may depend in any way on the matrices 7^^ and T". Note that in [2] this same approach 
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was followed for c** and c^, but not for the Dirac matrices 7^, which were replaced by 
ip^it). As a result the authors of ref. [2] obtained a supersymmetric sigma model, and the 
zero mode i/^o was useful to write the chiral anomaly in terms of p-forms in d dimensions. 
Since p-forms do not occur in trace anomalies, there is no advantage in replacing 7'^ by 
il)^[t) in our case. 

iv) If one defines the path integral by discretizing time into = 0, , . . . , , t = /3, 
the measure for q^{t) contains a factor which is due to integrating out the momenta 
p^{t) in the path integral. This factor is well-known from non-linear sigma models, and it is 
usually exponentiated, yielding a term with 5(0) which cancels the leading-loop divergences 
of the non-linear sigma model 0. In these factors were exponentiated by using scalar 
ghosts 

CgH = h{t)^Jg{q{t))c{t). (1.9) 
Below we shall instead use vector ghosts 

^9h = Ig^'iqit)) (6^(t)c^(t) + a^{t)a''{t)^ , (1.10) 

which give the same result, but which simplify the bookkeeping of loops because the qq 
part of the action is of the same form 

C^l9^''{q{t))(r{tW{t)]. (1.11) 



2' 

All divergences between the various non-ghost and ghost loops cancel for finite /3, which is 
rather remarkable as non-linear sigma models contain interactions with derivatives, con- 
firming the measure for q^{t). 

It is well-known that if one has intermediate points t^,t^, ■ ■ ■ ,tj^ between the initial 
and final times ti and tf, one needs (A^-fl) intermediate sets of momentum states, leading 



to (A^ + 1) factors y/del g [jl^. We shall exponentiate only A^ intermediate factors ^/g. 
As a result and vanish at the boundaries, and we can expand them into sin(7i7rt//3). 
The extra factor y/g appears then in the completeness relation as / dq^J g{q)\q>< q\ = 1, 
and leads to the transition amplitude 0, which is used to evolve wave functions 

i^{Qptf)= J dq,./^ <q^f,tM,U> i^{q>l,U) (1.12) 

v) In order to determine which Euclidean action corresponds to a given Hamiltonian 
(=regulator), there are in principle two ways to proceed. Either one starts with the Hamil- 
tonian operator formalism, discretizes time, introduces phase-space variables, integrates 
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out the momenta, and finally transforms from the discrete-time basis to, for example, a 
trigonometric basis. Or one bypasses this cumbersome but fundamental approach, and 
determines directly from the Schrodinger equation which Hamiltonian operator belongs to 
a given path integral. In the latter approach one starts thus with the Euclidean action 
and the path integral on a trigonometric basis, computes it to the required order in /3, and 
then finally finds the Hamiltonian to which it corresponds. We shall take a very general 
Euclidean action so that we can reach all the Hamiltonians we are interested in. The main 
problems with the former approach (the Hamiltonian operator approach) are the follow- 
ing, a) Given our regulator, there are no ordering ambiguities in the Hamiltonian itself; 
however, one has to evaluate the matrix elements <a;| exp{—eH)\p>, and here we should 
first expand the exponent, then move all x to the left and all p to the right, and then 
re-exponentiate the c- number result. For the harmonic oscillator this program is easily 
executed, but for non-linear sigma models it seems complicated, b) One has to compute 
the exact Jacobian for the transition from a discrete-time basis [the variables q{tj)] to 
a trigonometric basis [the coefficients qn in q{t) = '^qnSm{mTt/(3)]. Usually, one first 
guesses how the continuum action will look like, and then introduces a trigonometric basis 
whose measure one fixes such that it gives the correct result for the harmonic oscillator. 
Clearly, it would be better if no guesses or approximations were made, but we found this 
problem too difficult to solve. For further details, see appendix B. We shall thus follow 
Feynman's original method of deriving the Schrodinger equation, but not use the action 
which naively corresponds to the regulator, but rather a quite general action, which leaves 
room for extra terms. 

vi) The action we consider is of the form 



we must decide how to order the various terms in the expansion of exp(— 5"). If one uses a 
definition of the configuration space path integral where time is discretized, it is well-known 
that gauge invariance fixes an ambiguity with respect to the point at which we should take 
the gauge field. (The midpoint rule fixes the "Ito-ambiguity" , see ^^■) We shall expand 
all fields on a given trigonometric basis, and hence no ambiguity seems present. Yet, for a 




(1.13) 



The objects and V are matrices, and in order to define the path integral 




(1.14) 
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matrix-valued action, there is an ambiguity on how to order the terms in the expansion of 
exp{—S). We shall show that gauge covariance (not invariance) of our matrix- valued path 
integral selects time-ordering. As a mathematical exercise, we could have considered path 
integrals without time-ordering, but, as already said, this would violate gauge invariance; 
also, from a canonical operator point of view, time-ordering seems natural. 

In the following sections we shall describe these issues in detail. In sect. 2 we derive 
the regulators of the classical Weyl-invariant field theories, of which we shall compute the 
trace anomalies. In sect. 3 we define the path integral with matrix-valued connections 
and potentials. In sect. 4 we first expand the metric, connection and potential in terms of 
Riemann normal coordinates, and expand the quantum fluctuations of the coordinate and 
vector ghost fields in a Fourier series. We then evaluate the path integrals by making an 
expansion in terms of the parameter /3. We find the propagators and compute connected 
as well as disconnected Feynman diagrams on the world line. The result is a general 
formula for the trace anomalies in any dimension for abstract matrices and V in (4.22). 
Finally, in sect. 5, we apply this formula to the specific case of spin 0, |, and 1 in d = 2 
and d = 4, coupled to external gravitational and Yang-Mills fields. In sect. 6 we present 
our conclusions and comments. 



2. The regulators for classical Weyl-invariant theories 

In this section we obtain the regulators for the fields in the dimensions we are interested 
in. Before starting the discussion, we present our conventions. The Riemann curvature is 
defined by 

[v^,v,]v^^ = i^^,^v^^ (2.1) 

where is an arbitrary vector and the usual covariant derivative, which commutes 
with the metric. The Ricci tensor is given by = R^v^ <j and the scalar curvature by 

The action of a scalar field coupled to gravity and gauge fields reads as follows 



(2.2) 



where ^ = j^^r^j- The field (p transforms according to a given representation of the 
gauge group and (/) according to the conjugate one, so that (fxf) is a scalar under gauge 
transformations. The latter can be parametrized as follows 

A^^A^' = Ud^U-^ + UA^U-^; (f> ^ (f>' = Ucf), (2.3) 
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where it is understood that on the scalar field (/> the group element U acts in the corre- 
sponding representation (whose indices are suppressed). Under Weyl transformations 
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9fiu 9fiJ = ^g^jiu] A^' = A^; (j)' = fl 4 (p (2.4) 

the Ricci scalar transforms as 



R-^ R' = Q-^ 



R+id- 1) V^a^ log O + ^(rf - l){d - 2)g^^{d^ log Q){d, log O) 



(2.5) 



as a consequence the action ( p.2|) is invariant. The choice of TZ is not arbitrary if one wants 
to obtain a consistent anomaly A general method to identify such a regulator is to 
appeal to a Pauli-Villars regularization, which guarantees the consistency of the anomaly. 
For further details on this procedure, we refer directly to 0. For our purposes, we use 
the Fujikawa variables (p = g^cp, so that one is assured that there will be no gravitational 
anomalies, and we obtain the following expression for the regulator 

^0 = -Vi - ^R. (2.6) 

Here the gauge-covariant Laplacian acts in the representation which the field cp belongs 
to. If (p is real, then the anomaly is given by (|1.7| ), otherwise it is twice as large. 

The Weyl coupling of a complex (Dirac) spin | field ifj is described by the action 

Si = / dxe 'tpfipi 

(2.7) 

f = em^^^(d^ + u^ + A^, 



where e^"^ is a vielbein for the metric g^^, is its inverse, e = det e^"^ and 7"* are the 
SO{d) Dirac matrices. The spin connection = \oJ^rnnl^^: with 7"^"- = ^[7™', 7"'], takes 
values in the spinorial representation of the SO{d) Lie algebra and it is a function of the 
vielbein, while the gauge connection = A^T^ takes values in an arbitrary representation 
of an arbitrary Lie group. The Weyl transformations now read as follows 



Contrary to the scalar case, there is no need to add a curvature-dependent term to achieve 
invariance. To obtain the spin i regulator, we use again the fields ^jJ = g^i/j, and find 

= -y y = -vi,, -jR- ^F^.r" 

(2.9) 



A[^^ = -h lim Tr( ae"^^^ Y 
2 0^0 \ J 



Sf' = j d*x^ Tr 



For a real (Majorana) fermion one gets half this result. 

Finally we consider the spin 1 case, which is classically Weyl-invariant only in d = 4. 
We use the background- field method by splitting the gauge field into a background part 
and a quantum part A^f . Adding a gravitationally and background covariant gauge-fixing 
condition V'^A^f^^^ suitably weighted, one finds 

"''F.^F^^ + n (^v:^A(f") - ^n^ + Bv^^iv^c + [aI^''\c]) 

(2.10) 

where the subscript A indicates that the derivatives are background-gauge-covariant. We 
contract A^f^^^ with a vielbein field to obtain a tangent space index, and upon elimination 
of the auxiliary field 11, we find for the regulator of the spin 1 field: 

In addition we need the regulator for the ghosts. Both B and C are to be considered as 
real, and their regulator is 

Rbc = -V^. (2.12) 

These regulators act in the adjoint representation, thus the gauge connection as well as 
the Fmn term in eq. (|2.11 ) are absent for Abelian gauge fields. The trace anomaly is given 
by 



A^f^ = niim Tr( ae"^^-* ) - 2% lim Tr( ae 



-mi 



where we used 5 A 



(qu) 



(qu) 



SC = aC. and SB = 0, with A 



(qu) 



(2.13) 

9-^e^^Ai'-\ 



C = g^Cj and B = g^B. Perhaps we should note here that the gauge-fixing term in 
( ^.10| ) naively destroys the Weyl invariance. In fact, with the following natural definition 
of the Weyl transformations on the 11, B, C fields, 

n ^ n' = 

B B' = Q-^B (2.14) 
C ^ C" = C, 

the action is not invariant and there is no way of fixing the above transformation rules to 
achieve invariance. Under an infinitesimal Weyl transformation with Q = e'^ = l + a + ..., 
the action changes by 

n(<7^'^a.a)4^^) + Big^''d.a)iViC + [4'^^), C]) 



dSf^ = j d^x^ Tr 

= Sbrst J d'^Xy/g Tr 



(2.15) 



In the last line we have made use of the anticommuting BRST operator defined by 

SbrstC = — |C, C| 

2^ ^ (2.16) 

SbrstB = n 

SbrstTI = 0, 

which shows that the Weyl variation of the action is BRST-exact. For this reason it has 



zero expectation value between physical states [12 



3. Quantum mechanics with matrix-valued potentials in curved spaces 

We want to find now a path integral representation of the traces which yield the Weyl 
anomalies. The typical trace we have to compute has the following structure 

/ = lim Tr e~^^ 

(3.1) 

H = -{V''+An{d^ + A^) + S, 

where H acts on fields which are scalars under reparametrization and transform in a 
fixed, but otherwise unspecified, representation of the gauge group. The matrices 
and S transform as vectors and scalars under reparametrizations, respectively, but are 
furthermore quite general. Our strategy will be as follows. We first rewrite the trace in ( |3.1J ) 
using the operatorial formalism of quantum mechanics. Then we employ the path integral 
formulation of quantum mechanics to explicitly compute the trace. As is well-known, one 
can derive the path integral using the operatorial formulation of quantum mechanics and 
vice versa. As explained in the introduction, it is easier to start with a proper definition of 
the path integral on curved spaces, from which we derive unambiguously the Hamiltonian 
operator. 

The trace in (^]l|) can be taken over the Hilbert space generated by the following basis 
vectors: 

\x> (g)|z>= \x,i> (3.2) 

where the kets \i > correspond to a basis in the representation space of the gauge group 
and |x> are the usual kets of configuration space. The scalar product is defined by 

<x,i\yj>=^-^^^d,^ (3.3) 
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and the resolution of the identity takes the foUowing form 

i=j:fd.v-,\.,^><.,A (3.4) 

Any state | "0 > and operator G have the foUowing components in the above basis 

tpi{x) =<x, i\tp> 

(3.5) 

Gi^{x,y) =<x,i\G\y,j> 

In particular, the usual Hermitian position operator q^, Hermitian momentum operator 
and antiHermitian Lie algebra generator T" have the following matrix elements 

<x,i\r\y:j>=x^^-^^5^ 



X 



I- I ^ -1, ^ d 1. .S(x — y)^, 

<x,i\Pf,\y,3>= -^^9 ^{^)-Q^9^{x) (3-6) 



<x,tn,j>=T\ 



a ^(^ - y) 



\f9ix) 



One may check that is represented by —ihg~^d^g^ on wave functions iJi{x), and that 
it is Hermitian with respect to the inner product <(p\ip>= J dxg^(j)*ip. 

A gauge transformation is implemented by a unitary operator U and we have the 
following transformation properties for the wave function and the transition amplitude 



^|Ji{x) =<x, i\'i/j>-^<x,i\U\'ilj>= Ui-' {x)iljj{x) = i/^iix)' 
<x,i\e-^^\yJ>-^U,''{x) <x,k\e-^^\yJ>{U-')i^y), 



(3.7) 



where the Hamiltonian operator has the following form 

H = -rkip^. + A^Yg-^g^^^iip^ + i,)^-^ + S. (3.8) 
From (3.4) it follows that the trace in (|3.1| ) can be rewritten as 



J = lim ^ / dx^ <x,j\e ^^\xj> (3.9) 
J 

The gauge transformation properties in (^]^) for the transition amplitude as well as its 
scalar character under reparametrization will be the defining principles for our construction 
of the path integral representation, to which we now turn. 
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We are going to define the path integral representation for the transition amplitude, 
to be considered as a matrix with indices jf and jj, for a particle starting at point at 
the initial time ti and reaching point at the final time t / : 

{q';Jf,tM,J^,U) =<?y,i/|e-(*^-*^)^|Qf,j-> . (3.10) 

The path integral is defined as follows 

'•<litf)=qf 



{q^f,Jf:tM,j.,t,)= I ' '(^Q)Te-^M 

S[q] = j'j dt (^-g^.q^ct + + (3.11) 

{Vq)= n ^/9m)dq{t), 

ti<t<tf 

where T denotes time-ordering along the histories q^{t). The vector potential and the 
scalar potential V are matrices, so that the action itself is a matrix (hence the explicit 
indices ji and jf shown explicitly on the left-hand side of the first equation, but omitted 
in the right-hand side for notational convenience). The explicit relation between the scalar 



potential V in (|1.13| ) and (|3.11| ) and the scalar potential S appearing in (|3.1| ) will be given 
later. Note that the notation ( |...| ) means averaging with the functional integral under 
the boundary conditions given in the brackets, while < |...| > will denote the usual bras and 
kets of operatorial quantum mechanics, so that (|3.10| ) is the statement of the equivalence of 
the operatorial formalism and the functional integral formalism. To check the correctness 
of this definition of our path integral, we must check its transformation properties under 
reparametrizations and gauge variations. As far as reparametrizations are concerned, we 
can immediately see that the path integral is covariant (in fact it is a scalar), since it is 
built from manifestly covariant objects. Turning to the gauge symmetry, we can compute 
the variation of the transition amplitude under an infinitesimal gauge variation (U = 
with A infinitesimal) 

S{q1,jf,tf\q'^,ji,ti) = {q1,jf,tf\{A{qf) - A{qi))\q'^ , ji,ti) 

(3.12) 

= Aj/iqf){q>'f, k,tf\q^,j,, U) - {q'} J f,tf\q'^ , I, {q,). 

The first line is a consequence of the transformation properties of time-ordered exponen- 
tials. For the finite case it reads 

{q';,js,tAqt,h,U)' = U,,\qf){q^^,k,tf\q^^,tm~^V'{q^)^ (3-13) 
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which is the correct transformation property expected from the operatorial picture. 

We now turn to the path integral measure of ( |3.11| ). A key point introduced in 
is to exponentiate the non-trivial part of the measure using ghost fields, thereby leaving 
a translational-invariant measure necessary for performing the perturbative expansion. 
We present it here in a modified form, which helps in keeping track of the cancellation 
of potentially divergent terms. The latter are typically present in non-linear models of 
quantum mechanics (see ref. 0). Since we need to exponentiate the factor ^ present 
in the measure of (|3.11| ), we introduce anticommuting ghost fields h^^c^ and commuting 
ghost fields coupled to the metric g^^^. The effect of path integrating over h^^c^ will 
be to produce a factor g, while the path integration over will reduce it to thus 
recovering the correct measure. Therefore the covariant path integral which we are going 
to compute perturbatively in the next section reads as follows 

S\q, 6, c, a] = dt (^g^.M'^ct + ^''c^ + a'^a") + A^^cf -f (3.14) 

[Vq) = Yl dq{t): similarly for (Vb) (Vc) (Va). 

ti<t<tf 

4. Perturbation expansion 

We now present the perturbative evaluation of the covariant path integral defined 
in the previous section. Our strategy will be to make full use of the general coordinate 
invariance and pick the coordinate system best suited for carrying out the computation, 
namely the Riemann normal coordinates centered at the final point q'^. Before proceeding 
further, we introduce the dimensionless time variable r = — with /3 = tf — ti, and 
rescale the ghost fields as follows {h^^c^^a^) ^{h^^c^^a^)^ so that (3 will turn into a 
loop-counting parameter, nicely organizing the loop expansion. We do not keep track of a 
possible Jacobian for this rescaling, since the overall normalization will later be determined 
by other means. After such redefinitions the path integral takes the following form 

= r^''^~'\vq){Vh){Vc){Va) Te"^^!'^'^'^'"!, (4.1) 

J Q(ti )=Qi 



where 



.0 

S[q, 6, c, a] = / dr 
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(4.2) 



We now go over to Riemann normal coordinates by introducing functions (r, A) 
which are geodesies in A for fixed r, and which start at and end at q^ir): 

q>^{r, A = 0) = gy, q>^{r, A = 1) = q>^{T). (4.3) 

The normal coordinates z^^{t) for a point with coordinates q^ij) are then by definition 

(4.4) 



d 



z'^{T)^—q^{TA) 



A=0 



The action S[q{T, A = 1)] can be Taylor expanded in A around A = 0, and if the action is 
an Einstein scalar, one can replace ordinary A derivatives by covariant derivatives = 



n=0 ^ ^ 



(4.5) 



A=0 



Since the ^(r, A) arc geodesies, (^)^q(t. A) — 0, and in the expansion either no derivatives 
act on g^(T, A) (yielding the r- independent constant q^) or one derivative (yielding z^{t)). 
Another useful identity is 



D D 
dXd^ 



D D 



\T,X) = R^^^,z^q^z^, 



(4.6) 



where ^ = and z^{r, A) = ^q^{r, A). Of course, ^g^(r. A) = + 0(A). 

Using this expansion of the action into normal coordinates, one obtains 



)z^z''z^ 



+ i?-y5ap^/3(0)^^/^"^'' + ■ ■ (i^i" + 6^c" + a^^a") 

+ /3(^A^(0) + V.A^(O)^'^ + Ax,^{Q)z^z'' + ^«/?7m(0)^«^^^^ + 



(4.7) 



z^" 



+ (3^ (^F(O) + z'^V^F(O) + -z^z^V^s/^Vi^)) + 
where the following tensors have been defined 



1. 



1 



111 



(4.8) 
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Note the (3 dependence of the vector and scalar potentials produced by the introduction 
of the time variable r. It implies that they will start contributing at higher loops with 
respect to the "interactions" originating from the metric tensor. At this stage, we could 
start the perturbative quantum computation by splitting the field into a classical and 
a quantum piece 

= <i + <u- (4.9) 
The classical path z^^ is taken to be a solution of the classical equation of motion in the 
limit /9 — * 0, which is the geodesic equation. So z'^i{t) = —z^r, where z'^ are the Riemann 
normal coordinates of the initial point gf . It takes into account the boundary conditions 
prescribed to the path integral. The quantum fluctuations therefore vanish at r = —1,0 
and can be Fourier-expanded as follows 



^2^sin(7rnr). (4.10) 

n=l 

A similar expansion holds also for the ghost fields, which are zero at the boundary 

oo oo oo 

= ^ 6^ sin(7rnr); = ^ sin(7rnr); = ^ sin(7rnr). (4.11) 

n=l n=l n=l 

As we explained, these boundary conditions follow from the fact that the measure factor 
^/^ is not included in ( |3.11| ) for ti and tf. The evaluation of the transition amplitude in 
the two-loop approximation is enough to determine the quantum Hamiltonian associated 
with the path integral. Put differently, a two-loop computation is sufficient to verify that 
the transition amplitude satisfies a Schrodinger equation (diffusion equation) 

dt ^4_;L2) 
H = --(W + A^id^ + A,) ~^R + V. 

where in the first equation H acts on the indices and i. Such a computation was 
presented in and will not be reproduced here, since the introduction of the time ordering 
does not bring in any modification at this stage. The above Hamiltonian shows that the 
connection with the operatorial formulation presented in sect. 3 is achieved by setting 
the scalar potential V = + ^R. An additional outcome of the above procedure is the 
determination of the complete normalization of the path integral measure 



VqVhVcDa = -4 JJ JJ (7rm^)^(i2^d6^dc^(ia: 

m=l fJ,=l 
1 

(2^ 



(4.13) 



A = 

2 
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A consequence of this result is that all Gaussian integrals over the modes and the final 
integration over in (|1.12| ) are normalized to unity if one uses in ( [4.7| ) only the terms 
with gfj,u{0) and puts in ( |1.12|) the qi in ^/ g{qi) and il}{qi,ti) equal to qf. (The kinetic 
terms in the action are read off from (|4.7|) and are given by 

1 1 °° 

^(.) = 2^?M.(0)« + -^9,u{Q) ^^(Tr^n^^^^ + + a^^a^). (4.14) 

n=l 

The factors ^ due to integration over 2;^, and cancel, whereas the integration 

over z'^ produces a factor (27r/3)^((7)~i, which cancels against the {g)^ in (|1.12|) and A in 

(EH)-) 

Having established the equivalence (|3.10| ) for F = ^S" + we proceed to the task 
of computing the general trace in (|3.1J ) using the path integral representation 

/ = lim^ Tr e"^-^ = lim^ ^ j dxv^ (x^, z, 0|x^, z, (4.15) 



The evaluation of the right-hand side proceeds as follows. The classical path in (|4.9|) is 
now z'^i = 0, and we can drop the subscript in z^^ without possibility of confusion. From 
the quadratic piece in ( [1.14|) we recognize the following propagators 

{z^iT)z''iT')) = -(3g^''{0)AiT,T') 

{b^{r)c''{r')) = -/3^7^^(0)(2a2A(r,r')) (4.16) 
{a^{r)a''{T')) = -(3g^''{0){-d',Air,r')), 



where 

00 r r) 

A(r,r0 = 5: 

n=l 



sm Trnr sm tttit 



(4.17) 



Note that we can sum the series in ( [4.17| ) and obtain 

A(r, r') = r(r' + l)e(r - r') + r'(r + l)e(r' - r), (4.18) 

where the theta function is defined as usual by 

for T < t' 

e(r-r') = - for r = r' (4.19) 

1 for T > t' 
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Indeed, ^A(r,r') = 5{t - r'), while A(0,r') = A(r,0) = 0. Note that A(r, r') is 
continuous at r' | r and r' t and A(r, r) = r(r + 1). Also 

A(r, r') = •A(r, r') = r' + 6'(r - r') (4.20) 

(IT 

-^A(r, r') = A*(r, r') = r + ^(r' - r) (4.21) 
ar 

and also A* = 'A = ^ + r at r' = r. However, we must recall that the path integral 
is naturally regulated by truncating the mode expansion at a fixed mode M, so that the 
derivatives in ( [4.16| ) should be taken before the limit M — > oo. At fixed M, all propagators 
are well-defined functions, with finite limits at coinciding points t = t' . 

The interaction vertices are given by the remaining pieces of the action, S-^^ = S—S^^^ , 
with 5" given in ( [4.7|) and S^^,^ in ( [4.14|) . The final result for the basic trace is thus as follows. 
In d-dimensional Euclidean space, 

I^'^^ = lim Tr e"^^ = lim 3_ [ dx^(x)(0\Te~i^^^*\0) , (4.22) 

where 



(4.23) 



One must evaluate in d dimensions only all "vacuum polarization" graphs that are propor- 
tional to . They are obtained by expanding e'-"^ ^int\ and contracting all z^^b^^c^ 
and , using the propagators in ( [1.16| ) . The trace anomaly is then obtained by multiplying 
the result by the Weyl weights of the fields. 



5. Computation of the trace anomalies in d = 2 and d = 4 



To compute the trace anomalies in d = 2, /(^) in (|4.22 ), we need the (3 term in 
( [1.22|) . Since propagators go like (3 and vertices like ^ (or (3'^ for A, or (3 for V), this 
means that we must compute 2-loop diagrams in d = 2. Similarly for d = 2k we need the 
k + 1 loop contribution. Using the above propagators and vertices, we get for d = 2 only 
a contribution from the ^Ra^iui3 vertex in (|4.7| ). The corresponding diagrams have the 
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topology of the number 8, with none or one of the loops containing a ghost. The result is 
proportional to 



dr 



(•A* + A**)A- A'A* 



(5.1) 



The last term is clearly finite. Partially integrating the first term by using the fact that 
(•A* + A**) = 9t-(*A at r = r') immediately shows that ( |5.1| ) is finite. Also in higher loops 
finiteness is manifest thanks to using vector ghosts. The result is 



^(2) = try^^ 



(fx 
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(5.2) 



where the trace is over the finite dimensional representation space of the gauge group. 

The computation for obtaining the d = 4 result is more involved, but nevertheless 
straightforward. Let us first quote the result and then briefly comment on its derivation. 



!(4) 



tr. 



d X 
f2^ 



^ [R^^^pR'^^^P - R^.R^") + -^V^R 



720 



480 



1 

12 



^U^R 
2 V 24 



V 



+ —F F 
+ 48 



(5.3) 



The contribution due to one vertex —^S.^^ vanishes, since J^^dr •A(r, r) as well as 
J\dTA{T, r)*A(r, r) are zero. The contribution coming from two vertices —^S.^^ pro- 
duce the kinetic terms in -^F"^^. This comes about because / / drdr' {{z'^ z^){t) [z'^ zP){t')) 
is proportional to // drdT'lg""" g^'P /\{t, r')*A*(r, r') + g^Pg^'^'A^ir, r')*A(r, r')], and par- 
tially integrating 'A*, using that A vanishes at the boundary, gives the curl structure 
of the Maxwell action. The {DA)AA terms in the Yang-Mills action are a good test 
that our time-ordering prescription is required by gauge invariance. They are treated 
like the AV and AAV terms, see below. The V terms in (|5.3| ) come from the vertex 
-I3[V + ^z^z'^Vf.V^V]. First there is the term ^2^^[/^^]^ = 8^^^- ^he term with 
V'^V comes of course from {z^z'^)'V^'VuV. It is instructive to see how the covariantizations 



m 



Ld>^{d,V+[A„V]) 



j_ 

12 



A^^,d^V + [A^,V] 



arise, because here we see 



the crucial role of our time-ordering prescription. For example, from 

T(-^)2i2 /3(A^i^ + V,A^z^i^)rfr j'' ^(3^V + zPd,V)dT' 
zf'dT J zPdr' + /SdpVA^ J zPdr J z^'dr'. 



(5.4) 
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yields indeed —^[A'^.df^V]. Similarly 



one now has three 



Partially integrating and using z^{0) = z^{—l) 

the [d'^Afj^, V] is recovered. The most interesting term is A^, [A^, V] 

ways of ordering the time coordinates t^t'^t", and one indeed obtains the combination 
A^A^V - 2A^VA^' + VA^A^", and with the correct coefficient. 

We have been explicit in these calculations, because they show the importance of time 
ordering in order to obtain gauge-invariant results. 

Having computed the general traces (^T^) and ( [5.^ ) , we can specialize them to deduce 
the Weyl anomalies of the various fields discussed in sect. 2. We need to use V = ^R — ^^R 
for the spin regulator, and V = —\F^,^'~^^^^ and A^ A^ + uj^ for the spin i regulator, 
obtaining the following result for the trace anomalies of a real scalar and a complex spin 
^ field in two dimensions: 



In four dimensions we have instead 



{2n) 



:V9 



instead 

aR^^^pR'^^''^ + bR^^R^"" + cR^ + dV'^R + eF^^F^"" 



a 



(5.5) 



(5.6) 



with the various fields contributing as follows 
spin a = b 

spin I (Dirac spinor) a = 
spin 1 (Abelian) a = 



1 

720 



b=^ 
^ 180 



720 " 90 



C 
C 






d = 


1 

720 


e = 


1 

48 


1 

288 


d = 


1 

120 


e = 


1 
6 


5 
144 


d = 


1 

40 


e = 





5 
144 


d = 


1 

40 


e = 


11 

24 



(5.7) 



spin 1 (non- Abelian) a = —j^ ^ ~ M 
where in addition we have used for the ghost regulator V = ^R and for the gauge regulator 
V = ^R5m^ — \Rm^ — |-Fm" and A^ — > + W/^, with uj^ acting in the vector representa- 
tion. These results agree with the literature . For completeness, we recall that eq. ( |5T6| ) 
can also be written in the following form 

d'^x 
(2^ 

where JF = Ra^upR"'^^^ — 2R^i,R^^ -\- ^R^ is the square of the Weyl conformal tensor and 
G = RctiupR"'^"'^ - AR^^R^""" + R^ is the topological Euler density. The coefficients a, /3, 7 
are given by 



aT + l3g + -fV'^R + e(tr 



(5.^ 



1 



a 

P 
7 



480 



(A^o + 6A^i + 12Ni 



1 



1440 
1 



(A^o + llA^i + 62A^i 



(5.9) 



720 



(A^o + 6A^i - 18A^: 
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where A'^o, A^i and A^i are the number of real spin 0, complex (Dirac) spin i, and real spin 
1 fields, respectively. This form of presenting the trace anomaly clearly shows that the 
coefficients in eq. ( |5.6|) are not all independent 0. This is a consequence of the Wess- 
Zumino consistency conditions, which guarantee the absence of a single term in the 
anomaly fl^. Moreover, the coefficient 7 in ( ^.81) can be modified at will (in particular 
can be set to zero) by adding to the effective action a finite local counterterm proportional 
to R^. For spin 1, only the sum of the contributions of the non-ghost sector and the ghost 
sector satisfies the consistency conditions. The gauge-fixed action is not Weyl-invariant, 
but its Weyl variation is BRST exact and has vanishing expectation value ||12||. 



6. Conclusions 

In this article we have computed trace anomalies in quantum field theory by using 
path integrals in quantum mechanics. Our work extends the analysis of Alvarez-Gaume 
and Witten of chiral anomalies, but because trace anomalies are not of a topological nature, 
certain fundamental and interesting aspects of path integrals had to be taken into account, 
which could be neglected in the case of chiral anomalies. As already indicated in |T^, this 



makes the evaluation of trace anomalies more complicated than that of chiral anomalies. Of 
interest is the method itself, rather than the results. The quantum mechanical approach is 
an excellent tool for computing anomalies. It avoids the painful Baker-Campbell-Hausdorf 
expansion of a direct application of Fujikawa's method. It is clearly related to the DeWitt's 
heat-kernel approach [^, but it avoids the bilocal tensor calculus of the latter. 

We have given a general formula, eq. (|4.22|) , from which trace anomalies in any 



dimension can be obtained. We then computed the trace anomalies in d = 2 and d = 4 for 
spin 0, ^ and 1 at one loop with external gravitational and Yang-Mills fields. Our results 
agree, of course, with the literature 0] . In the literature also the trace anomalies for d = 6 
and d = 8 have been obtained (by Gilkey [|16[ and Avramidi [H^, respectively, who both 
used the heat-kernel approach). It would be interesting to find an explicit closed formula 
for the trace anomalies in any d, similar to the result of 0] based on the "Dirac genus" for 



the chiral anomalies. Perhaps our eq. ( 4.22 ) could be a starting point. 

In the case of the chiral anomaly, one needs in d dimensions d powers of the fermionic 
zero modes i/jq in order that the Berezin integration is non-vanishing. This introduces a 



factor P^, which cancels the factor /3 2 in eq. ( |4.13| ). For the trace anomaly, no i/'o are 



present, and hence one must make a loop expansion on the world line. 
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In the transition from phase space to configuration space path integrals, a measure 
^Jg is produced, famihar from non-hnear sigma models. We have treated this measure in 
a way that is different from the usual 5(0) approach 0. Namely, rather than add extra 
vertices to the action, which are proportional to In det g times 5(0) and which are supposed 
to cancel leading divergences, we exponentiated this ^ to a perfectly regular term in the 
action with new ghost fields |]^, similar to the familiar Faddeev- Popov ghost action. It 
was then clear, diagram by diagram, how all the divergences (not only the leading ones) 
cancel in all loops for finite /3. An example of such a cancellation was given in ( |5.1| ). Of 
course, we need to exponentiate the factor ^ in order to be left with a translationally 
invariant measure from which propagators can be defined. Historically, these factors 5(0) 
first occurred in the proof of Matthews 's theorem fl^, which states that the Hamiltonian 
and the Lagrangian canonical quantization are equivalent if one adds in the latter extra 



vertices proportional to 5(0), and in the analysis of Lee and Yang of massive vector 
bosons. We believe that our approach is simpler and clearer. 

In the canonical formalism, Heisenberg fields satisfy the field equations, and fields 
must be expanded in a complete set of solutions. In the path integral, fields do not need to 
satisfy any field equations, and one is free to expand in any basis. We found it convenient to 
decompose the coordinates q^ify into a background piece, which satisfies the field equation 
(a geodesic), and fiuctuations. We expanded all fiuctuations in the trigonometric basis 
sin(7rmr) with m = 1, 2, 3, . . .. This raises the question whether our results depend on the 
choice of basis. Since an ambiguity is known to exist if one uses time discretization, and is 
fixed by imposing gauge invariance, we expect our results to be basis-independent as long 
as they are Einstein- and gauge-invariant. 

We end with a comment on Feynman diagrams, which provide the oldest method 
of computing anomalies and which have been used with success for chiral anomalies 0. 
The simplest diagram which one might consider to yield a trace anomaly is the two-point 
function in d dimensions, with candidate trace anomaly (V^)^~^-R. In d = 2 it is a genuine 
anomaly, but for 2(i> 2 it can be removed by a counterterm with two factors K (as observed 
by Adam Schwimmer (private communication) the most general two-point function which 
is both transversal and traceless is proportional to \{d — \){Qp,pQva + ^^.a^up) — dfiudpa^ 
where 9^i, = r]^^^ — k^ki,k~'^, but in d = 2 this expression identically vanishes). The 
triangle diagram cannot reproduce all h? terms in (i = 4, because the h? terms of the 
Euler invariant vanish in all dimensions. Thus, in (i = 4, one would need to compute 
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the 4-point function, which is tedious to evaluate. Moreover, in higher d there are more 
invariants, which constitute the trace anomaly, and hence further diagrams would have to 
be computed. (For the chiral anomaly, only one diagram needs to be calculated in a given 
dimension.) We conclude that the quantum mechanical approach is a much simpler way 
to compute the trace anomalies. 
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Appendix A. 

In the computation of loops on ttie world line, we have used the explicit form of the 
propagator function A(r, r'). As a result, all loop integrals became trivial, typically powers 
of r and r' from —1 to 0. However, if one is interested in truncated results one needs to 
use the representation of A(r, r') in terms of sine functions, and then the loop integrals for 
d = 4 become (when the truncation is removed at the end of the calculation) proportional 
to the following double sums 

^ (m + nYmn 180' {m-\-n)m?n 72' ^—^ {m-\-nYm? 120 

n,m=l n,m=l n,m=l 

2 

The sums Yl^=i given by the polylogarithms Li2p-2(1); for example Li2(l) = 

Li4(l) = 1^, Li6(l) = ^ and Li8(l) = gf^, see [|^. Since Zl^^i ^"^ = T"' factor tt"^ 



in ( |A.1D is plausible, and one easily finds the coefficients and from a numerical 

approximation. However, since (to our surprise) the results for these double sums are not 
found in the standard references [^, we give an analytical derivation, which can also be 
generalized to more complicated sums. 

First, consider /(p, r) = X]^m=i ("^ + '^)~^"^~^'^~'^- -^y decomposing the sum- 
mands for fixed n into partial sums with m-dependent singularities and n-dependent coef- 
ficients, one derives the following identities 

/(I, 2, 1) = /(I, 0, 3) + /(O, 2, 2) - /(0, 1, 3) 
/(2, 1, 1) = -/(2, 0, 2) - /(I, 0, 3) + /(0, 1, 3) 

(A.2) 

/(3, 1, 0) = -/(3, 0, 1) - /(2, 0, 2) - /(I, 0, 3) + /(0, 1, 3) 
/(I, 3, 0) = -/(I, 0, 3) + /(O, 3, 1) - /(O, 2, 2) + /(0, 1, 3). 

(From /(2, 2, 0) one obtains the same result as from /(I, 3, 0).) From these relations, one 
finds /(I, 2, 1), /(2, 1, 1) + /(2, 2, 0) and 2/(3, 1, 0) + /(2, 2, 0), which agree with dOl) . 
To obtain also /(2, 1, 1), A. Martin gave us the following derivation 

^-^ [ra + nYran ^ — ^ /n x /n y m n 

n,m=l n,m=l 

= r^[-ln(l-,)P= r^ln^(l-,) r^ = - r^ln^(l-:.)lnx. 

Jo ^ Jo y Jo y Jy X Jo ^ 

(A.3) 
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fl — x)"^ — 1 1 , r 1 T.2e_ c 

, but tor in x we use 



to cancel the - singularity. 



Next we replace ln(l — x) by 

One then obtains a series of beta functions B{x,y) = dt t^~^{l — t)y~^ — Y(x+y) 
being the gamma function: 



r(x)r(y) 



7(2,1,1) = -^ 



S(2e, 1 + 2e) - S(e, 1 + 2e) - 2S(2e, 1 + e) + 2S(e, 1 + e) + ^ - - 



• (A.4) 



Next we use r(2e) = ^'^^^^'^^ and obtain 
/(2,1,1) = 



1 

2? 



2e 

r2(l + 2e) \ /r2(l + e) r(l + e)r(l + 2e) 



r(l + 4e) 



-1+4 



r(l + 2e) 



r(l + 3e) 



(A.5) 



Writing r(l + z) as exp(z'i/'(l) + ^z'^ijj'{l) + . . .) where tp{z) = P^^ (digamma) 

function, given also by iIj'{z) = Xl^o + ^)~^^ ^^^^ V'(l) ^ind all even deriva- 

tives of i/j{z) at z = 1 cancel. The 2n + 1 derivatives of iIj{z) at z = 1 are proportional to 
C(2n), which are known: = J2n=i^~^ = ^'"(1) =^En=i^~^ = f^- Substitu- 

tion of these expansions into 7(2, 1, 1) then reveals that all e poles cancel, and the finite 
result is /(2, 1, 1) = 

We thank A. Schellekens for discussions, and expecially A. Martin for a derivation of 
7(2,1,1). 



Appendix B. 

In this appendix we comment on some problems that are one encountered in the 
operator approach to anomalies. We start from operators Ai, . . . , An depending on 
canonical variables a = ^^^^^ and = ^^^^ , satisfying [d, d^] = 1. [If the Hamiltonian 
is i7 = + ^mw^Q^, then q = {mu)^Q and p = {mu)~^P, so that dPdQ = dpdq 

and H = (d^d + Using the so-called holomorphic representation [^, [^, the trace 

of the product of the Ai can be written as a multiple integral (where dwdw = idpdq, etc.) 

^ , - - . f dwdw %-r dvjdvj ( _ _ \ . , x a i- \ 

Tr(Ai,...,Aiv) = / / Qy.\>\-ww - ^ViVi\Ax[w,vx) . . .AN{yN-\,w). 

j=i ^ i=i ' 

(B.l) 

The kernels A{z,z) are equal to (exp Z2;)i^(^, 2;), where K{z,z) = Xlm n>o 

n>o^ -^"ind". Note that this last expression has to be normal ordered with 
respect to d^ and d. 
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Consider now exp i-PH) = (exp(-ei7)) , where A^e = (3. Let exp(— ei?(a^^, a)) = 
: exp(—eh{a^, a)) : where h will in general depend on e, then the trace can be written as 

N-l 



_ fiTT f dwdw -i-r dvidvn , ^, 



' ^ _ (B.2) 

-L = -w{w - Vi) - vi{vi -V2) + ...- vn-i{vn-i - w) 

- t[h{w, vi) + h{vi, V2) + ... h{vN-i,w)]. 

By identifying w = vn and lu = vq, this result can be written in a suggestive form as a 
path integral with periodic boundary conditions (PBC) 

but the exact meaning (exact even for finite N) is given by (|B.2|) . 

Given the operator H (the regulators in our case), it is very difficult to find h{z, z) in 
closed form. One has to expand exp(— eff), commute all to the left and d to the right, 
replace d^ by z and d by z, and finally re-exponentiate. For the harmonic oscillator this 
has been done but not, to our knowledge, for non-linear sigma models. For the harmonic 
oscillator with H = [a) d + \)uj one obtains 



exp(— e/i(2, z)) = exp 



--ew — exp(l — e '^^)zz 



(B.4) 



To obtain the configuration space equivalent of ( p.2| ), we must integrate out the N 
momentum variables p,pi, . . . ,p7v-i- In the continuum limit (|B.3|) we find from the term 
Jq drv^, upon re-discretization, a sum Yl^=o ^(Qj ~ ^Pj)iQj + ^Pi)^- Owing to PBC, this 
can be rewritten as a sum of j-^iOj + Pj)^ — WjQj^ and if h{p,q) is quadratic in p [so 
h{p,q) = \g^-' {q)piPj + V{q)] one obtains the Euclidean action L = dT{^gijcfq^ + V) 
together with the measure ^ discussed in sect. 1. 

However, if one starts from the exact result in (p.2|) , the computation becomes con- 
siderably more complicated. Consider first the harmonic oscillator. One obtains then in 
the exponent an expression of the form — | ^(q'j — ipj)Mjk{qk + ipk)^ — \Ntuj, where the 
non-symmetric matrix M has entries along the diagonal, and —e~'^'^ just above the 
diagonal and in the lower-left position. The integration over the variables pj yields 

(27r)*(detMs)-^ / TT ^ exp f- Jg(Ms - M^M^ 1m^)q') (B.5) 
J PBC jj^ zn \ z J 
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and further integration over the q variables yields the expected result detM = 1 — 
exp(— ewA^). Together with the factor (exp— ^ew)^, one finds (2sinh^)~^, which is 
indeed the partition function for the harmonic oscillator. To see detM coming out, note 
that 

deiMs det(Ms - MaMs^Ma) = deiMs{Ms - MA)Mg\Ms + Ma), (B.6) 

where Ms = ^{M + M^) and Ma = ^{M — M^). For a system more general than 
the harmonic oscillator, one would have to compute detMg in order to obtain the exact 
corresponding configuration space path integral, and this is in general difficult. (Even for 



the harmonic oscillator, we have not found detMg. In [jT^ the determinant is given of a 
symmetric matrix, which is equal to our Ms except that the lower left entry and upper 
right entry are zero.) 

For the chiral anomaly, one must evaluate Tr exp(i^F) exp(— where F = c^c is 
the fermion number operator. The integral kernel for Ai = exp{i9F) can be found in 
closed form, because the variables c and c are anticommuting: Ai(c, c) = exp(e*^cc). The 
path integral corresponding to (|B.3|) is obtained by replacing the first factor exp{—eH) by 
exp(i6'-F), and one finds 

N-l 



^WF^-PH ^ /■ « TT ^MLexp(-L) 
J 2ni -LJ. 2ni ^ 



-L = - eA.-i) + . . . + - + ee^'(ei + e-^'O 

-e[/i(6,6) + ---/i(?"iv-i,0]- 

We have taken care to define ( |B.1D such that it also holds for anticommuting variables, as 
in ( [13.71 ). By identifying = C and = — e~*^^ we obtain a path integral as in ( [tj.3[ ), but 



now with antiperiodic boundary conditions ABC for ^ = 0, and PBC for 9 = n. Hence, 
for the trace anomaly we would need fermions with ABC, see sect. 1. 

Finally, we consider the ghosts c*^ and q for internal symmetry, discussed in sect. 1. 
We have found the projection operator Pi onto the one-particle states c*^|0 > in closed 
form. It is given by 

Pi =: exp(c^^Ci) : -1. (B.8) 

Written out. Pi = c^^q + ^c'-^c^^CjCi + . . ., and it is clear that Pi|0>= and Pic*^|0> = 
c*^|0 >. One can check that Pi annihilates the states with 2, ...,M particles (where 
j = 1, . . . , M). It follows that the analysis is very similar to that of the chiral anomaly: 
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one must only replace the integral kernel exp(e*^cc) by exp(cc) — 1. However, we have not 
been able to interpret this result in terms of suitable boundary conditions. 

Our conclusion is that the operator approach, which starts from the quantum Hamil- 
tonian, is much more complicated than the Feynman approach, which produces the Hamil- 
tonian at the end. For this reason we have followed in the main text the Feynman approach 
(but we have not made the error of omitting quantum fluctuations from the path integral 
p5[1 ). The approach of the authors in was based on an operator approach, but they 
made various approximations: they assumed that h = H, and used ( [B.3| ) instead of ( |B.2| ), 
to obtain the Euclidean action. For the chiral anomalies, their approximate path integral 
yields the correct result because chiral anomalies depend on very few details of the path 
integral (they are topological). For the trace anomaly these approximations would be 
incorrect, as we discussed in the introduction. 
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